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ABSTRACT 


The  data  handling  and  analysis  of  recorded  d.ita  from  . 
group  of  Independent  buoys  deployed  In  the  vicinity  of  the 
Chesapeake  Light  Tower  Is  described.  Included  Is  a 
description  of  the  digitization-quality  checking  process* 
the  analysis  technique  used  and  a discussion  on  th<»  results 

SUMMARY 

The  David  W.  Taylor  Naval  Ship  Research  and  Development 
Center  (DTNSRDC)  has  developed  software  to  digitize  and  ana'yze 
U S Coast  Guard  Exposed  Water  Buoy  data  tapes.  A total  of  48 
analog  magnetic  tapes  from  the  initial  buoy  deployment  (Dec  .’2  - 
Jui  73)  have  been  logged*  quality  checked*  digitized  and  analyzed 
Each  tape  contains  frequency  multiplexed  motion  and  tension  signals 
from  three  moored  buoys  as  well  as  wave  height*  time  code,  voice 
and  calibration  Information.  The  analysis  Is  based  on  data  from 
3$  runs  selected  for  quality  and  duration  and  Is  representative 
of  sea  states  2 and  3.  The  analysis  results  Include  both  tabular 
and  graphical  output,  such  as  power  spectral  densities,  response 
amplitude  operators,  cross  spectrums  and  statistical  Items, 
resulting  directly  from  the  data  analysis  computer  programs  and 
graphical  displays  derived  from  these,  such  as  run  to  run 
comparisons,  peak  amplitude  resets  and  averaged  results.  Also 
included  are  X-Y  plots  of  some  of  the  raw  analog  data. 

Section  1 of  this  report  describes  the  data  logging,  quality 
checking  and  digitization  procedures  as  well  as  some  data  acquisition 
background  information.  Section  2 Includes  a description  of  the 
analysis  techniques  used.  Section  3 discusses  the  results  of 
applying  this  analysis  to  the  subject  data  and  Appendix  A is  the 
foreword  for  the  computer  generated  analysis  results  which,  because 
of  the  large  volume,  are  bound  separately 


INTRODUCTION 


Buoys  have  long  been  used  as  a marine  navigational  aid.  The 
effectiveness  of  a given  buoy  design  Is  a function  of  Its  cost,  ease 
of  handling,  and  Its  usefulness  to  the  mariner,  e.g.,  the  detectability 
and  recognition  of  Its  signal.  Evaluation  of  alternative  designs  and 
comparison  with  the  standard  8'  x 26'  buoy  requires  quantitative 
Information  on  these  factors.  Since  mooring  costs  and  detectability 
are  both  greatly  affected  by  the  motion  characteristics  of  the  buoy, 
motion  measurements  are  necessary  to  provide  quantitative  information 
for  design  evaluation. 

Consequently,  a buoy  motion  measurement  system  was  designed  and 
Installed  to  provide  buoy  motion  and  driving  function  data  for  evaluating 
navigational  buoy  design  concepts  In  the  operational  environment.  A 
future  goal  of  this  program  will  be  the  Implementation  of  a dynamic 
simulation  of  a moored  buoy  In  a regular  and  irregular  seaway.  Tabulated 
and  plotted  results  from  the  data  analysis  will  be  used  for  comparison 
with  simulation  generated  data.  The  primary  objective  of  this  overall 
effort  Is  to  aid  In  selecting  the  most  effective  designs  for  new 
buoy-bssed  navigation  systems  used  by  all  vessels  in  coastal  waterways. 
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SECTION  1 - BACKGROUND 


1-1  Data  Acquisition 

The  U.  S.  Coast  Guard  (USCG)  Buoy  Motion  Measurement  System  was 
developed  by  the  David  W.  Taylor  Naval  Ship  Research  and  Development 
Center  (DTNSRDC).  Wave  height  at  the  light  tower  and  roll  angle,  pitch 
angle,  heave  displacement  and  mooring  line  tension  from  up  to  seven  buoys 
located  near  the  tower  are  measured  and  recorded  on  command.  Although  the 
system  can  accommodate  up  to  seven  buoys,  only  buoys  number  3}  4 and  7 
were  active  during  the  collection  of  -ata  for  this  report.  Buoys  numbers 
4 and  7 are  the  standard  8'  x 26'  configuration  and  buoy  number  3 Is  the 
OSPREY  design  (see  Figure  1).  Note  that  buoy  number  3 is  moored  by  a 
single  chain  from  Its  central  mooring  eye,  while  the  standard  buoys  are 
bridle  moored.  The  weight  In  service  of  the  OSPREY  is  approximately 
7700  ibs  and  of  the  standard  approximately  5000  lbs. 

A block  diagram  of  the  entire  USCG  Buoy  Motion  Measurement  System  Is 
shown  in  Figure  2.  The  buoy  Instrumentation  Is  shown  on  the  left  side  of 
the  diagram.  Buoy  instrument  packages  are  Installed  at  the  roll  and  pitch 
center  of  rotation  and  are  powered  by  batteries  located  In  the  buoy.  Buoy 
instrument  packages  are  commanded  to  turn-on  and  calibrate  by  the  command 
unit  located  at  the  ground  station.  Buoy  roll,  pitch,  and  heave  are 
measured  by  a gyroscope  and  a vertically  stabilized  accelerometer  located 
in  the  buoy  package.  Mooring  cable  tension  Is  measured  by  means  of  an 
external  tension  dynamometer  which  Is  located  In  the  mooring  cable  and 
electrically  connected  Into  the  package.  The  conditioned  data  from  the 
four  transducers  is  multiplexed  (mixed  together)  and  transmitted  to  the 
data  receiver  located  at  the  ground  station. 

The  ground  station  equipment  Is  shown  In  the  center  of  the  block 
diagram.  The  tape  recorder  Is  controlled  by  the  command  system, 
Multiplexed  data  from  each  of  the  seven  buoys  Is  recorded  separately 
on  one  of  the  fourteen  tape  channels.  In  addition  to  buoy  data,  wave 
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Figure  2 - System  Block  Diagram 


information  from  a Bailor  Wave  Profile  System  located  on  the  outside  of  the 
tower,  time  code,  a reference  frequency*  an  "OK  data"  signal,  and  voice 
Information  are  also  recorded  on  the  tape  recorder. 

The  data  tapes  are  transported  to  the  data  analysis  center,  represented 
on  the  right  of  Figure  2,  where  the  information  for  each  buoy  is 
de-multiplexed,  quality  checked,  digitized  and  analyzed  on  a digital  computer. 
The  time  code  translator  Is  utilized  to  locate  times  and  dates  on  the  data 
tapes. 

1*2  Analog  Tapes 

DTNSRDC  Is  furnished  magnetic  tapes  containing  data  recorded  In  analog 
form  on  a Honeywell  Model  5600  Instrumentation  tape  recorder.  The  recorder 
uses  a standard  one  Inch  magnetic  tape,  TRIG  standard  record  heads  and 
the  data  was  recorded  at  1 7/8  Inches  per  second.  Data  and  information 
from  the  light  tower  (ground  station)  is  recorded  on  Channels  1-7  ("A" 

Inputs).  The  Channels  5 and  7 are  reserved  for  water  current  not  presently 
recorded.  Channels  8 through  14  ("8"  Inputs)  are  reserved  for  the 
subcarrier  multiplexed  data  from  seven  buoys.  Each  multiplexed  channel 
contained  8 subcarrier  measurements.  The  complete  data  set  on  the  analog 
tape  is  described  In  Table  1.  Table  2 describes  the  "A"  Inputs  In  detail 
and  Table  3 describes  the  "8"  inputs. 

The  analog  tapes  contain  the  wave  height  sensor  data,  data  from  the 
three  buoys,  and  Identification  and  control  information.  The  time  code 
channel  Is  used  to  Identify  Individual  data  runs.  The  OK  data  channel 
serves  to  indicate  the  beginning  and  ending  of  data  runs  and  to  Identify 
calibration  signals.  Each  of  the  experimental  data  channels,  the  wave 
height  and  the  four  buoy  parameters  are  recorded  at  two  different 
sensitivities.  The  appropriate  sensitivity  is  selected  during  the  data 
reduction  process. 

The  analog  data  is  reproduced  at  a tape  speed  of  30  inches  per  second, 
sixteen  times  real  time.  An  eight  channel  demodulator  and  filter  are  used 
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to  de-multiplex  the  recorded  subcarrier  buoy  data  Into  8 separate  analog 
data  channels  acceptable  for  Input  Into  the  digital  processing  system. 

1-3  Data  Logging  and  Quality  Checking 

Analog  tapes  were  received  from  the  Coast  Guard  at  a rate  of  about 
5 tapes  per  week  during  the  period  between  June  and  September  1973.  Input 
logging  of  the  analog  tapes  is  a three  step  process.  Figure  3 shows  a 
typical  Coast  Guard  Data  Recording  Log  which  is  received  with  each  analog 
tape.  The  analog  tape  Is  assigned  a Tape  Number  and  an  DTNSRDC  Tape  Log 
Is  prepared  (see  Figure  4).  The  starting  time  codes  listed  on  the  Coast 
Guard  Log  are  transposed  to  the  DTNSRDC  log  and  run  numbers  are  assigned 
corresponding  to  tape  numbers  and  time  code.  During  the  Quality  Checking 
and  Digitizing  Process,  the  DTNSRDC  log  is  completed  to  Include  the  date 
digitized,  the  number  of  good  records  and  quality  remarks.  Note  that  each 
analog  tape  will  result  In  one  DTNSRDC  log  per  buoy  (presently  three). 
Figure  5 shows  an  analog  Tape  Status  Summary  which  was  established  to 
provide  the  Coast  Guard  with  continuing  feedback  relative  to  DTNSRDC 
progress  during  the  digitization  and  analysis  process. 

Upon  receipt  of  a group  (usually  5 to  10)  of  analog  tapes,  a strip 
chart  record  of  the  most  recent  tape  Is  generated  and  delivered,  with 
comments,  to  DTNSRDC  Code  296  for  determination  and  analysis  of  problems 
in  the  data  acquisition  equipment.  This  early  notification  to  the  Coast 
Guard  point  of  contact  as  to  the  condition  of  the  analog  tapes  permits 
timely  Coast  Guard  preventive  action  for  future  tapes. 

Initial  data  quality  checks  on  the  analog  tapes  are  performed  by 
visual  observation  on  a strip  chart  recorder  In  parallel  with  the 
digitization  process  (dotted  line  on  Flgifre  6).  The  strip  chart  record 
of  the  next  run  Is  examined  during  the  digitization  process  of  the  present 
run.  Parallel  digitization  and  analog  quality  checking  with  operator 
manual  override  Is  more  efficient  than  performance  of  these  operations 
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TABLE  1 - TAPE  RECGRDER  DATA  CHANNELS 


Tape 

Record 

Data 

hannel 

Electronics 

Recorded 

VA 

Voice 

Remarks* 

1 

Direct 

Time  Code* 

2 

FM 

OK  DATA* 

3 

Direct 

Reference 

4 

FM 

25  FT  Wave  Height 

5 

FM 

SPARE 

6 

FM 

5 FT  Wave  Height* 

7 

FM 

SPARE 

8 

Oirect 

Buoy  #1 

9 

Direct 

Buoy  02 

10 

Direct 

Buoy  #3 

11 

Direct 

Buoy  #4 

12 

Direct 

Buoy  #5 

13 

Oirect 

Buoy  #6 

14 

Direct 

Buoy  07 

* Wave  Height  Only  (WHO)  Recording  System 


8 


TABLE  2 - "A"  INPUTS  FOR  ANALOG  TAPE 


(Voice)  Time,  date,  run  number,  visual  observations  of 
wave  height  (P-P),  wave  period  (sec),  and  general  remarks. 

(Direct)  IRIG  B time  code  (days,  hours,  minutes,  seconds). 

Log  should  Indicate  the  approximate  time  code  for  the 
beginning  of  each  run. 

(FM)  "OK  DATA"  signal.  Indicates  when  calibrating  and 
when  recording  data. 

(Direct)  6.25  KHZ  reference  signal  used  for  tape  speed 
compensation  in  the  discriminators  on  playback.  Playback 
Is  at  30  IPS  (16  x 1 7/8  IPS)  which  results  In  a compensation 
frequency  of  100  KHZ. 

(FM)  50  ft  full  scale  (+  25  ft)  wave  height  signal. 

(FM)  spare  channel  reserved  for  ocean  current. 

(FM)  10  ft  full  scale  {+  5 ft)  wave  height  signal,  clips 
signals  greater  than  + 5 ft. 

(FM)  spare  channel  reserved  for  ocean  current. 


>^j^-  v ■ , ^ i 


TABLE  3 - HB"  INPUTS  FOR  ANALOG  TAPE 


CH8-14-  (Direct)  multiplexed  buoy  date  on  IRIG  telemetry  channels  1-2. 


IRIG 

CH 

MEASURAND 

-CAL  (OV) 

0 CAL 
(2.5V) 

♦ CAL 
(5.0 V) 

1 

100°  F.S.  (+50°)  ROLL 

-50°  STBO* 

0° 

♦50°  PORT* 

2 

20°  F.S.  (+10°)  ROLL 

-10°  STBO* 

0° 

+10°  PORT* 

3 

100°  F.S.  (+50°)  PITCH 

-Su-  AFT* 

0* 

+50°  FWD* 

4 

20”  F.S.  (+10°)  PITCH 

-10*  AFT  ' 

0° 

♦10°  FWD* 

5 

1 G F.S.  (+.56)  VERT.  ACCEL. 

-.5  G OWN 

0 G 

+.5  G UP 

6 

.2  G F.S.  (+  1G)  VERT.  ACCEL. 

-.1  G DWN 

0 G 

♦.1  G UP 

7 

50  KIP  F.S.  TENSION 

0 KIP 

25  KIP 

50  KIP 

8 

10  KIP  F.S.  TENSION 

0 KIP 

5 KIP 

10  KIP 

* Direction  of  rotation 
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accusations 
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USC6  >nalo9  Tape  Status  Summary 


Hgur*  6 - Analog  T«p«  Procosslng 


sequentially.  The  following  criteria  Is  used  for  deciding  which  runs 
contain  good  Information  for  a sufficient  duration  to  warrant  further 
analysis 

. calibration  signals 
. gyros  fully  erected 
buoys  turned  on 
steady  zero  lines 
. excessive  noise  and  spikes 
. duration 

. discontinuities  In  the  data 

If  the  calibrating  signals  are  bad,  the  digitization  process  has  the 
capability  of  providing  artificial  calibrations.  When  artificial  calibrations 
are  inserted  it  is  noted  under  remarks  on  the  OTNSRDC  log  sheet  Also 
listed  on  the  log  sheet  is  the  operator  s estimate  of  the  sea  state  in  1 FT 
peak-to-peak  Increments  of  wave  height  observed  from  the  wave  height  signal 
on  the  strip  chart  recorder.  If  the  data  observed  on  the  strip  chart  does 
not  warrant  digitization  then  a note  is  made  In  the  DTNSRDC  log.  Normally, 
bad  data  Is  referred  to  In-one  of  two  categories;  (1)  NO  G000  - data  did 
not  meet  examination  criteria  or  (2)  NO  DATA  - no  data  on  the  tape  for 
the  indicated  time-code 

A final  quality  check  Is  made  at  tne  end  of  the  digitizing  process 
when  a read  and  print  of  the  digital  tape  is  made  and  the  operator  inspects 
the  printout  Assuming  the  tape  passes  this  final  inspection  it  is 
transferred  to  the  analys’s  computer  for  processing  The  entire  analog 
tape  processing  procedure  just  described  is  illustrated  >n  Figure  6 

1-4  Digitizetlon 

A block  diagram  of  the  digitization  system  for  the  processing  of  the 
analog  tapes  Is  shown  in  Figure  7 The  analog  tape,  which  was  recorded  at 
1 7/8  inches  per  second,  is  reproduced  at  30  IPS  (16  times  real  time)  on 
an  Ampex  FR-1800  tape  recorder.  An  eight  channel  demodulator  is  utilized 
to  de-multiplex  the  individual  data  from  each  buoy  to  etght  separate  ana  leg 
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Tower  FM  Data 


Figure  7 - Block  Diagram  of  the  Processing  Equipment  for  the 
Digitization  of  Analog  Tapes 


data  signals  (2  sensitivity  levels  for  each  of  the  four  measurands).  The 
demodulator  Is  equipped  with  a tape  speed  compensation  unit  A 100  KHZ 
reference  signal  recorded  at  6.25  KHZ  Is  utilized  for  this  compensation  The 
time  code  generator  is  used  to  properly  locate  and  identify  each  run. 

The  de-mu lti pi exed  buoy  signals  along  with  the  remainder  of  the  tape 
recorder  signals  are  fed  through  6 Pole  Butterworth,  Low-Pass  filters  to 
remove  high  frequency  flunctuatlons  and  noise.  The  filter  half-power  point 
is  set  at  16  HZ  which  reflects  a real  time  frequency  of  1 HZ  (6  28  rad/sec) 
The  filter  outputs  are  simultaneously  directed  to  the  A/D  converter  and  the 
strip  recorder.  The  A/D  converter  converts  the  analog  signals  to 
digital  format  and  feeds  them  into  the  computer.  The  computer  is  commanded 
via  the  teletype.  The  digitized  data  Is  stored  on  magnetic  tape  and  ID 
arrays  are  outputted  to  the  printer. 

The  first  data  recorded  during  a data  run  Is  a series  of  calibration 
levels  These  levels  are  Identified  by  three  negative  excursions  of  the 
OK  data  signal  - the  first  Indicating  negative  cal,  the  second  zero  cal 
and  the  third  positive  cal.  The  experimental  data,  consisting  of  the  two 
wave  height  channels  and  the  eight  signals  from  the  buoy  cf  interest,  are 
digitized  at  discrete  time  intervals-.  The  portion  of  the  actual  data  to 
be  digitized  is  defined  by  a positive  excursion  of  the  OK  DATA  signal  The 
digitized  data  and  calibrations  are  recorded  on  an  intermediate  (scratch} 
tape.  The  scratch  tape  Is  next  re-read  back  into  the  computer.  The  program 
selects  the  data  range  to  be  used  for  each  variable  based  on  the  number  of 
times  the  lower  range  exceeds  full  scale.  The  three  calibration  levels 
for  the  selected  data  channels  are  averaged  to  provide  three  single 
calibration  values  for  each  variable  The  plus  and  minus  calibration  levels 
are  summed  to  provide  a single  peak-to-peak  calibration  value  for  each 
variable-  The  averaged  zero  calibration  level  is  used  to  convert  for  zero 
offset  or  drift  in  the  data  signals.  The  five  data  signals  (variables) 
selected  and  corrected  along  with  the  ID  arrays  (final  peak-to-peak 


calibration  values,  run  numbers,  no.  of  records,  record  no.,  buoy  no., 
and  Julian  date)  are  outputted  and  stored  on  the  final  digitized  tape. 

The  following  procedure  is  followed  In  digitizing  one  run  for  one  buoy 
from  an  analog  tape 

(1)  load  analog  tape  to  start  of  run  (time  code  and  voice  is 
used  to  assist  m this  operation) 

(?)  based  on  prior  inspection  of  the  run,  select  or  not  select 
artificial  calibration  (switch  on  computer) 

(3)  load  scratch  and  data  tapes  on  digital  tape  recorders 

(4)  initialize  computer,  make  proper  logical  unit  assignments 
and  input  via  teletype  - run  number,  sampling  rate,  buoy 
no.  and  Julian  date 

(5)  start  tape  recorder 

(6)  watch  OK  data  line  on  strip  chart  recorder,  when  OK  data 
goes  negative  for  the  first  calibration  - start  the  computer 

(7)  after  analog  tape  has  completed  run,  let  recorder  continue 
while  digitization  process  Is  being  completed  and  quality 
check  the  data  for  the  next  run 

(8)  log  in  no.  of  records,  no.  of  good  records,  remarks  and  sea 
state  in  ft  peak-to-peak 

1-5  Digitization  Status 

To  date,  a total  of  48  analog  tapes  containing  1857  runs  from  the 
Initial  buoy  deployment  (Dec  72  - Jul  73)  have  beer  received,  logged  and 
quality  checked.  Of  these  1857  runs,  a total  of  1126  were  digitized  after 
quality  checking  The  731  runs  not  digitized  were  eliminated  for  the 
following  reasons  - no  data  on  the  tape,  gyros  not  fully  erected,  buoys  not 
turned  on,  unsteady  zero  lines,  excessive  noise  and  spikes,  Insufficient 
duration  or  discontinuities  in  the  data.  Further  inspection  of  the  data 
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mealed  that  the  wave  height  signal  on  an  additional  524  runs  suffered  from 
a periodic  noise  Interference  signal  and  thus  these  runs  were  of  no  further 
value  for  analysis.  Therefore*  602  runs  remained  as  candidates  for  analysis. 
Of  these,  109  runs  contained  good  data  for  all  3 buoys  over  a significant 
Interval  of  time  (l.e.,  greater  than  10  minutes).  Table  4 shows  a breakdown 
of  these  runs  with  respect  to  visually  observed  wave  height  and  number  of 
good  records.  The  25  runs  containing  at  least  100  records  (greater  than 
14  minutes)  were  selected  for  further  analysis. 


TABLE  4 - BREAKDOWN  OF  109  GOOD  RUNS 


# GOOD  RECORDS 

r 

VISUALLY  OBSERVED  WAVE  HEIGHT 
2'  3'  4'  5' 

6' 

136-140 

mm 

mm 

«•  m 

mm 

131-135 

mm 

mm 

1 

1 

mm 

m m 

126-130 

m m 

2 

m m 

mm 

121-125 

— 

1 

m m 

mm 

m m 

m m 

116-120 

2 

1 

1 

m m 

mm 

mm 

111-115 

4 

m m 

1 

m m 

mm 

106-110 

3 

2 

1 

mm 

mm 

m m 

101-105 

1 

2 

2 

m m 

— 

mm 

> 100  TOTALS 

10 

8 

6 

1 

m m 

— 

96-100 

mm 

12 

6 

m m 

mm 

91-95 

1 

4 

2 

mm 

mm 

mm 

86-90 

3 

4 

4 

— 

mm 

— m 

81-85 

m m 

7 

• m 

m m 

m m 

mm 

76-80 

2 

8 

m m 

mm 

mm 

71-75 

2 

2 

1 

1 

— 

m m 

> 70  TOTALS 

18 

45 

19 

2 

m m 
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SECTION  2 - ANALYSIS  TECHNIQUES 


Two  analysis  techniques  are  used  in  the  data  analysis  computer  module 
of  the  OTNSROC  Buoy  Motion  Measurement  System.  These  are  Statistical  Analysis 
and  Spectral  Analysis  the  details  of  which  are  discussed  in  this  Section. 

2-1  Statistical  Analysis 

The  wave  height,  buoy  vertical  acceleration,  pitch  and  roll  angles, 
and  the  mooring  cable  tension  were  analyzed  as  follows: 

Let  X.j,  1 * 1 , 2 ... , M,  be  the  time  sequence  digitized  values  of 
one  of  the  5 variables  and  M Is  the  total  number  of  values.  The  following 
statistics  can  then  be  computed: 

M 

Average  x * £ x 'M 

1*1  1 


Mean  Square 


£ x.  /M 
1*1  1 


Variance  C.  £ x.^/M  - x2 

° i«i 1 

Standard  deviation  o * JCT 

o 

Probability  density  * frequency  of  occurrence/M 

Probability  distribution  at  b * E frequency  of  occurrence  at  b 

n n 

for  all  b„  > b 
n - 


In  most  instances,  the  full  scale  of  a variable  is  divided  Into  50  sections. 
The  frequency  of  occurrence  of  each  section  bn  Is  tallied  in  order  to  compute 
the  probability  density. 
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For  the  sea  wave  the  significant,  height  Is  defined  as 


* 4 C 
o 


where  CQ  Is  the  variance  of  the  wave  amplitude. 

2-2  Correlations  and  Power  Spectral  Densities 

The  cross-covariance  of  two  variables,  x and  y,  for  k number  of 
lags  h (sample  periods)  Is  computed  as 


1 M-k 

cxyk  * rrr  ^ vi+k  • * J 


The  auto-covariance  of  x Is  the  cross-covariance  of  x and  x Itself: 


1 M"k  2 
Cxyk  ' fmr  x1x1+k  -x 

Auto  or  cross-covariance  Is  the  unbiased  (or  zero  biased)  auto  or  cross- 
correlation. 

The  power  spectral  density  Is  computed  In  the  following  manner: 
Assume  that  the  auto-covariance  (or  cross-covariance  In  the  case  of 
cross-spectral  density)  Is  becoming  less  and  less  correlated  as  k 
Increases  and  Is  practically  zero  when  k approaches  N.  N is  an 
Integer  less  than  H.  (As  a rule  of  thumb,  N should  be  less  than  M/10). 
Consider  then  the  auto-correlatlon  Is  periodic  with  period  T « (2N  + 1) 
h from  - (N  <*  1/2)  h.  The  Fourier  series  of  this  periodic  function  at 
kh,  k * 0,  1,  2,  ... , Is 
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c,  <kh)  ■ »0/2  ♦ L,  t *p  cm  fSfftrb.  * P s,n  liiftNH  J 


V P*»  v P 


with  the  coefficients 


a * 2 C (t)  cos  p w t dt  p ■ 0»  1,  2, 
*p  * S xx  ' o 

» T/9 


b * t / (t)  sin  p w t dt 

P ‘-T/2  ** 


and  the  fundamental  frequency  Is 

“o  * m ?rnr 

Both  a and  b„  art  computed  In  th»  way  described  by  Goertzel  In  Ralston 
and  wolf's  Mathematical  Methods  for  Digital  Computers,  pp.  25B-262. 


The  spectral  density  of  x Is  defined  as 

ell  1 7 C„  (t)  e"lut  dt 

S„  M • 5,  ..  « 


Under  the  assumption  that  C„  Is  practically  zero  for  large  k the 
spectral  density  «t  w * pwQ  1* 

T/2 

J*  (P»0)  ■ 1 Jm  C„  (t)  C cos  pujp  Sin  P-0  t ] dt 

. T C ap  - 1 bp  ] 
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The  auto-covariance  Is  an  even  function  and  bp  * 0.  For  cross-spectrum 
sxy  th*  expression  Is  the  same  If  It  Is  computed  from  the 

cross-covariance  Cxy.  However,  Is  not  an  even  function  and  bp  Is 
not  zero,  although  the  assumption  that  Cxy  Is  practically  zero  for  k close 
to  and  beyond  + N has  to  be  applied  In  the  same  way. 

Among  the  buoy  variables  (pitch  and  roll  angles,  vertical  acceleration 
and  cable  tension  on  the  buoy)  the  frequency  response  transfer  function 
between  two  variables  x and  y Is 

H <*•>  * ffr^- 

The  cross-spectral  density  S^x  Is  a complex  function  of  u.  It  can  be 
expressed  In  polar  form  with  amplitude  Am  and  phase  a. 

H (J“>  * rTuT " 8 

The  spectral  density  obtained  In  this  way  Is  discrete,  only  the 
values  at  puQ,  p * 0,  1,  2,  ...,  N,  are  defined.  The  fundamental 
frequency  uQ,  therefore,  determines  the  resolution  of  the  spectrum. 

Suppose  uQ  Is  taken  to  be  .1  rad/sec  or  .1/2*  * .0159  Hz.  The 
fundamental  period  Is  then  1.0159  - 62.89  seconds.  Since  the  sample 
period  of  the  data  points,  or  the  lag,  Is  .16  seconds,  the  fundamental 
period  Is  then  corresponding  to  62.89/  16  * 382  lags.  Therefore,  to 
have  a .1  rad/sec  resolution  of  the  spectrum,  the  auto  or  cross-covariance 
must  be  computed  over  382  lags.  Since  the  auto-covariance  Is  an  even 
function.  It  Is  computed  for  k « 0,  1,  2,  ....  196  only  (H  ■ 196).  For 
cross-covariance,  we  let  k range  from  -196  to  +196. 

Finer  resolution  can  be  obtained  with  wave  lags,  but  the  computation 
of  Cxxk  or  Cxyk  very  consum^n9*  W*  consider  the  fundamental 
frequency  «Q  * .1  rad/sec  1$  adequate  and  use  N » 196. 
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Due  to  the  assumption  of  zero  correlation  for  k greater  than  N and 
the  finite  period*  the  spectral  values  at  0 and  .1  rad/sec  are  not 
accurate  in  actual  practice.  To  smooth  out  the  spectrum*  a three  point 
"hannlng"  process  Is  also  applied  to  the  spectral  density  values.  Because 
of  this*  the  point  at  .2  rad/sec  Is  not  accurate  either.  Therefore,  the 
first  couple  of  points  (w  < .2  rad/sec)  In  the  spectrum  should  not  be  used. 

2-3  Response  Amplitude  Operators 

For  the  sea  wave  with  wave  frequency  w (rad/sec)  and  length  X,  the 
deep  water  relation  Is: 


The  mean  wave  slope  is  defined  as  2;/X,  where  c Is  the  wave  amplitude. 
The  wave  slope  Is  non-dimensional  and  Its  spectral  density 


J M“)  * Try  S (w) 
x*  4 4 


Is  used  tc  compute  the  non-dimensional  RAO  (response  amplitude  operator) 
of  the  roll  and  pitch  angles.  The  S^(u)  Is  the  sea  wave  amplitude 
spectral  density  with  amplitude  In  feet. 

2 

If  Sr  (w)  Is  the  roll  angle  spectral  density  In  degree  - seconds, 
the  RAO  for  roll  is  defined  as 


57. 3u 


1.754  • Sr(h)* 

ur  s>) 

•*  c 


The  same  applies  to  the  pitch  angle,  replacing  Sf  (w)  by  the  pitch 
spectral  density  (w). 

The  RAO  of  buoy  vertical  acceleration  Is  normalized  to  sea  acceleration 
and  Is 


where  S (u)  Is  the  buoy  vertical  acceleration  spectral  density  with 
» 2 
acceleration  In  ft/sec  . 
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SECTION  3 - ANALYSIS  RESULTS 


This  Section  describes  the  results  generated  directly  from  the  data 
analysis  computer  module  of  the  DTNSRDC  Buoy  fetion  Measurement  System, 
as  well  as,  those  results  derived  Indirectly  from  these.  Interpretation 
of  these  results  and  conclusions  based  on  them  are  also  discussed. 

Twenty-five  runs  were  used  In  the  analysis  (see  Table  S).  Note  that 
for  each  run,  over  100  records  were  processed  with  respect  to  each  of  the 
three  buoys.  The  output  generated  directly  by  the  computer  is  bound 
separately,  and  Is  provided  with  Appendix  A.  It  consists  of  a set  of 
spectral  analysis  results  and  a set  of  statistical  results.  The  75 
pairs  (for  each  run,  3 pa1rs  corresponding  to  the  3 buoys)  of  these  sets 
are  presented  In  the  order  In  which  they  appear  In  Table  5. 

At  the  beginning  of  each  set  of  spectral  analysis  results,  there  are 
two  tables  of  data  upon  which  the  plots  that  comprise  the  remainder  of 
the  set  are  based.  The  top  of  the  first  table  contains  two  rows  of 
summary  data  that  are  self-explanatory.  The  remainder  of  this  table 
consists  of  the  power  spectral  density  (PSD)  for  the  5 parameters: 
wave  height  (indicated  as  sea),  roll  angle,  pitch  angle,  heave  displacement 
(indicated  as  accel),  and  tension  and  the  response  amplitude  operator 
(RAO)  for  the  latter  four  parameters.  Both  the  PSO's  and  their  RAO's 
correspond  to  frequencies  between  .1  and  5 radians/sec  Incremented  by 
1 radian/sec.  The  second  table  contains  the  following  5 Items 
(tabulated  for  the  frequencies  Just  mentioned)  representing  the  cross 
spectrum  of  heave  displacement  and  pitch  angle  and  the  cross  spectrum 
of  heave  displacement  and  tension:  In  phase  term.  In  quadrature  term. 

In  magnitude.  In  phase  angle,  and  transfer  function,  Reference  Is  made 
to  Section  2 for  the  mathematics  related  to  the  computation  of  the 
PSD's,  RAO's  and  cross  spectra. 

Each  set  of  statistical  results  Is  headed  by  run  summary  data  almost 
identical  to  that  heading  the  spectral  analysis  results.  Plots  of  the 
probability  density  of  wave  height,  heave  acceleration,  total  angle. 
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roll  angle,  pitch  angle  and  tension  are  each  preceded  by  related  statistics 
such  as  maximum  value*  minimum  value,  average  value,  root  mean  square 
value,  etc.  Again,  reference  Is  made  to  Section  2 tor  the  mathematics 
related  to  these  statistical  computations. 

The  Indirectly  derived  results  can  be  thought  of  as  falling  into 
three  categories  - run  to  run  comparisons,  peak  amplitude  results,  and 
averaged  results.  The  run  to  run  comparisons  (Figures  8*31)  Mere 
generated  by  plotting  results  from  3 consecutive  runs  on  the  same  set 
of  axes.  Two  groups  of  such  results  were  plotted,  one  corresponding  to 
sea  state  2 (runs  36.14,  37.05,  37.11)  and  one  corresponding  to  sea  state 
3 (runs  36,12,  38.16,  36.13).  In  each  group  there  are  12  plots:  Wave 
Height  vs  Frequency,  RAO  of  Roll  Angle  vs  Frequency  (one  plot  for  each 
of  the  3 buoys),  RAO  of  Pitch  Angle  vs  Frequency  (one  plot  for  each  of 
the  3 buoys),  RAO  of  Heave  Displacement  vs  Frequency  (one  plot  for  each 
of  the  3 buoys),  and  Phase  Angle  vs  Frequency  of  the  Cross  Spectrum  of 
Heave  Displacement  and  Pitch  Angle  (for  buoys  3 and  4).  The  Have  Height 
vs  Frequency  plot  is  not  related  to  any  one  buoy  since  the  wave  height 
measurement  Is  independent  of  the  buoy.  It  is  Interesting  to  note  that 

on  each  of  the  RAO  plots,  the  values  for  the  3 runs  are  quite  close. 

Further,  the  roll  angle  plots  corresponding  to  buoys  4 and  7 (the 
Standard  Buoys)  are  almost  identical,  whereas  this  similarity  Is  not 
as  strongly  exhibited  with  respect  to  pitch  angle  and  heave  displacement. 

In  these  latter  two  cases,  the  RAO  values  corresponding  to  buoy  7 appear 
to  be  greater  than  those  for  buoy  4 over  large  segments  of  the  frequency 
range  (e.g.,  from  2.0  to  3 8 rad/sec  for  heave  displacement,  sea  state  2). 

The  peak  amplitude  results  (Figures  32-49)  consist  of  the  peak  RAO 

of  roll  angle,  pitch  angle  and  heave  displacement  for  each  of  the  3 buoys 

plotted  as  a function  of  wave  height  standard  deviation  and  then  density 
at  wave  peak.  A single  plot  corresponds  to  a single  buoy  and  one  of  the 
three  buoy  motion  parameters  versus  either  wave  height  standard  deviation 
or  density  at  wave  peak.  Each  plot  contains  data  from  three  sea  states, 
the  points  corresponding  to  standard  deviations  greater  than  1.30  (or 
wave  peaks  greater  than  1.60)  are  from  sea  state  4 runs,  the  point 
corresponding  to  standard  deviations  less  than  .75  (or  wave  peaks  less 
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than  .25)  are  from  sea  state  2 runs  and  all  the  In  between  points 
are  from  sea  state  3 runs. 

The  averaged  results  consist  of  plots  of  the  means  (Indicated  by  x) 
and  90%  confidence  limits  for  the  RAO  of  roll  angle*  pitch  angle  and  heave 
displacement.  Two  groups  of  results  are  plotted*  one  corresponding  to  sea 
state  2 (Figures  50*58)  and  one  corresponding  to  sea  state  3 (Figures 
59-67).  A single  plot  corresponds  to  one  of  the  two  sea  states*  one 
of  the  three  buoys,  and  one  of  the  three  buoy  motion  parameters.  For  a 
given  frequency*  the  mean  Is  computed  by  averaging  the  RAO  values 
corresponding  to  that  frequency  over  the  runs  corresponding  to  the 
given  sea  state.  For  sea  state  2 (see  Table  5)  the  runs  averaged  over 
36.14  to  37.15  Inclusive  and  for  sea  state  3,  38.07  to  39.07  Inclusive. 

The  901  confidence  limits  are  given  by 

s 

5 ♦ (t)  n 

where  x * sample  mean  (l.e.,  the  average  over  the  appropriate  runs) 
s « sampel  of  standard  deviation 
n « number  of  runs  used  to  compute  x (or  s) 
t * Student's  t-dlstrlbutlon  value  corresponding  to  n*1 
degrees  of  freedom  and  901  confidence  level 

One  other  category  of  results  Is  Included  here  In  an  attempt  to  give 
the  reader  a better  visual  representation  of  the  buoy  motion  * analog  x-y 
plots  of  roll  and  pitch  angles.  The  analog  tape  was  played  back  at  1 7/8 
Ips  through  real  time  discriminators  and  then  Into  an  x*y  recorder  adjusted 
for  appropriate  scaling  - 1°  per  minor  division  and  10°  per  major  division. 
Pitch  angle  Is  plotted  on  the  900/270°  axis  and  roll  angle  on  the  0°/180° 
axis.  The  graphs  therefore  display  50°  full  scale  for  roll  angle  anJ  35° 
full  scale  for  pitch  angle.  Three  runs  were  plotted  - 26.10,  37.04  and 
39.04.  For  run  36.10,  time  segments  of  1.25  minutes  (Figures  68-70),  2.50 
minutes  (Figures  71-73),  5 minutes  (Figures  74-76),  10  minutes  (Figures 
77-79)  and  20  minutes  (Figures  80-82)  were  plotted.  For  the  other  two  runs, 
only  10  minute  time  segments  were  plotted  (Figures  83-88).  In  each  case, 
all  3 buoys  are  shown. 
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Figure  8 - Wave  Height  versus  Wave  Frequency  In  Sea  State  Two 
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Figure  10  - RAO  of  Roll  Angle  versus  Wave  Frequency  for  Buoy  4 
(8  x 26  - 08)  In  Sea  State  Two 
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Figure  12  - RAO  of  Pitch  Angle  versus  Wove  Frequency 
for  Buoy  3 (OSPREY)  In  Sea  State  Two 
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Figure  20  - Wave  Height  versus  Wave  Frequency  In  Sea  State  Three 


BUOY  3 
SEA  STATE  3 


WAVE  FREQUENCY  (RADIANS/SECOND) 


3T9MV  110M  JO  OV* 


44 


WAVE  FREQUENCY  (RADIANS/SECOND) 


Figure  25  - RAO  of  Pitch  Angle  versus  Wave  Frequency  for 
Buoy  4 ( 8 x 26  - 08)  in  Sea  State  Three 
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Figure  26  - RAO  of  Pitch  Angle  versus  Have  Frequency  for 
Buoy  7 ( 8 x 26  - 05)  in  Sea  State  Three 
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Phase  Angle  versus  Wave  Frequency  of  the  Cross  Spectrum  of  Heave 
Displacement  and  Pitch  Angle  for  Buoy  3 (OSPREY)  In  Sea  State  Three 
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Figure  32  - Peak  RAO  of  Roll  Angle  versus  Wave  Height 
Standard  Deviation  for  Buoy  3 (OSPREY) 


Peak  RAO  of  Roll  Angle  versus  Wave  Height 
Standard  Deviation  for  Buoy  7 (8  x 26  - 05) 


Figure  35  - Peak  RAO  of  Roll  Angle  versus  Spectral  Density 
at  Wave  Peak  for  Buoy  3 (OSPREY) 


5»«r73T. 


Figure  37  - Peak  RAO  of  Roll  Angle  versus  Spectral  Density 
at  Wave  Peak  for  Buoy  7 (8  x 25  - 05) 


Figure  38  - Peak  RAO  of  Pitch  Angle  versus  Wave  Height 
Standard  Deviation  for  Buoy  3 (OSPKtf) 


Figure  41  - Peek  RAO  of  Pitch  Angle  versus  Spectral  Density 
at  Wave  Peak  for  Buoy  3 (OSPREY) 


Figure  43  - Peek  RAO  of  Pitch  Angle  versus  Spectrel  Density 
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Figure  44  - Peak  RAO  of  Heave  Displacement  versus  Have  Height 
Standard  Deviation  for  Buoy  3 (OSPREY) 
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Figure  48  - Peek  RAO  of  Heave  01spUce*«nt  versus  Spectral 
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Figure  49  - Peek  RAO  of  Heave  Displacement  versus  Spectral 
Density  at  Have  Peak  for  Buoy  7 (8  x 26  - 05) 
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WAVE  FREQUENCY  (RADIANS/SECOND 
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Figure  62  - Neon  Velue  of  RAO  of  Pitch  Anglo  versus  Wove 

Frequency  for  Buoy  3 (OSPREY)  In  See  Stete  Three 
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FREQUENCY  (RADIANS/SECOND) 


Figure  U - Anelog  X-i  Plot  of  Roll  «nd  Pltc*  Angles 
for  Buoy  3 (OSPREY) 


of  Roll  end  Pitch  Angles 


Figure  70  - Aneloc  X-Y  Plot  of  Roll  end  Pitch 
for  Buoy  7 (8  * 26  - 05) 
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Figure  74  - Analog  X-Y  Plot  of  Roll  and  Pitch  Angles 
for  Buoy  3 (0SPR£Y) 


X-Y  Plot  of  Roll  and  Pitch  Ang)** 


Figure  77  - Analog  X-Y  Plot  of  Roll  and  Pitch  Angles 
for  Buoy  3 (OSPREY) 


Figure  78  - Analog  X-Y  Plot  of  Roll  and  Pitch  Angle* 
for  Buoy  ♦ (8  x 26  - 08) 


Figure  79  - Analog  X-Y  Plot  of  Roll  and 
for  Buoy  7 (8  x 28  - OS) 
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Figure  80  - Analog  X-Y  Plot  of  Roll  and  Pitch  Angles 
for  Buoy  3 (OSPREY) 
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Figure  85  - Analog  X-Y  Plot  of  Roll  a ni  Pitch  Angles 
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Figure  87  - Analog  X-Y  Plot  of  Poll  ami  Pitch  Angles 
for  Buoy  4 (8  x 26  - 08) 
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APPENDIX  A - COMPUTER  GENERATED  ANALYSIS  RESULTS 


The  David  W.  Taylor  Naval  Ship  Research  and  Development  Center 
(DTNSRDC)  has  developed  software  to  digitize  and  analyze  U.  S.  Coast 
Guard  Exposed  Water  Buoy  data  tapes.  The  analysis  Includes  data  from  a 
group  of  Independent  buoys  deployed  In  the  vicinity  of  the  Chesapeake 
LT  Tower  during  the  Initial  buoy  deployment  (Dec  72-Jul  73).  Each 
tape  contains  frequency  multiplexed  motion  and  tension  signals  from  up 
to  seven  buoys  as  well  as  wave  height*  time  code*  voice  and  calibration 
Information. 

Although  the  system  can  accommodate  up  to  seven  buoys,  only  buoys 
number  3,  4 and  7 were  active  during  the  collection  of  data  for  this 
deployment.  Buoy  numbers  4 and  7 are  the  standard  8’  x 26’  configuration 
and  buoy  number  3 is  the  OSPREY  design.  Note  that  buoy  number  3 is 
moored  by  a single  chain  from  its  central  mooring  eye,  while  the 
standard  buoys  are  bridle  moored.  The  weight  In  service  of  the  OSPREY 
is  approximately  7700  lbs  and  of  the  standard  approximately  5000  lbs. 

The  analysis  Is  based  on  data  from  25  runs  selected  for  quality 
and  duration  and  Is  representative  of  sea  states  2 and  3.  These  25 
runs  contain  at  least  14  minutes  of  continuous  data  for  buoys  3,  4 
and  7.  The  data  Is  divided  Into  75  separate  files  (3  buoys  x 25  runs) 
and  Is  tabulated  In  order  of  decreasing  variance  of  wave  height  In  Table 
A1 . The  analysis  results  contain  both  tabular  and  graphical  output  and 
Include  75  pairs  of  statistical  and  spectral  analysis  results. 

In  addition  to  mooring  cable  tension,  buoy  motions  Incl'de  roll, 
pitch  and  heave  acceleration  all  referenced  to  the  vertical.  Wave 
height  was  measured  at  the  light  tower  and  recorded  along  w ,h  the 
buoy  motions  and  reference  Information. 
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In  addition  to  basic  statistics  (max,  Min,  average,  RMS,  MEANSQ, 
variance  and  sigma)  the  statistical  segment  also  includes  plots  of 
probability  amplitude  densities  and  a probability  density  and  distribution 
of  total  angle  Spectral  Analysis  results  include  tables  and  plots  of 
power  spectral  density,  Response  Amplitude  Operators,  Cross  Spectrums 
including  magnitude  and  phase  and  Transfer  functions.  Some  abbreviations 
used  are  as  follows: 

SPS  - samples  per  second 
FS  - full  scale  sensitivity 
CAL  - calibration 
PSD  - power  spectral  density 
RAO  - response  amplitude  operator 
SEA  - wave  height 
CSURFACE  - wave  height 

ACCEL  - heave  acceleration  (statistical  segment) 

ACCEL  - heave  displacement  (spectral  segment) 

A detailed  description  of  the  digitization  - quality  checking  process, 
the  analysis  technique  used  and  a discussion  on  the  results  are  included 
In  OTNSRDC  Report  No  SP0-657-01. 
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| TABLE  AT  - TWENTY-FIVE  RUNS  USED  IN  THE  ANALYSIS  - RUN  #'s  LISTED  IN 

I ORDER  OF  DECREASING  VARIANCE  OF  WAVE  HEIGHTS 


Sea 


* 

State 

Run  # 

Buoy  #3 

Buoy  #4 

Buoy  #7 

f 

4 

55.11 

5.3189 

• 

5.3616 

l 

4 

38.05 

1.8193 

1.8198 

1.8131 

r» 

u 

5* 

4 

38.06 

1.7595 

1,7226 

1.7091 

f 

5 

3 

38.07 

1.5499 

2.0648 

1.5009 

3 

36.10 

1.5406 

1.5170 

1.5444 

L* 

f 

3 

36.11 

1.3740 

1.8293 

1.9727 

£ 

3 

36.12 

1.3545 

1.3347 

1.2821 

[ 

3 

38.16 

1.1884 

1.1948 

1.1897 

r 

t 

3 

36.13 

.8267 

.8464 

.8350 

f 

\ 

3 

39.04 

.7977 

.7941 

.7981 

if 

3 

37.04 

.6580 

.6466 

.6541 

3 

41.06 

.6559 

.6465 

.6458 

3 

39.07 

.6450 

.6402 

.6656 

f 

2 

36.14 

.5249 

.5358 

.5300 

< 

2 

37.05 

.5095 

.5003 

.5129 

' 

* 

2 

37.11 

.4751 

.4676 

.4805 

V 

? 

2 

37.02 

.4508 

.4373 

.4502 

f 

f 

2 

36.15 

.4427 

.4447 

.4527 

c 

: 

* 2 

37.12 

.4083 

.3994 

.4073 

■ 

2 

36.06 

.3716 

.3418  ' 

.3500 

V 

L 

2 

38.01 

.3657 

.3554 

.3645 

[ 

2 

36.16 

.3562 

.3505 

.3597 

2 

37.14 

.3592 

.3504 

.3570 

I 
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OTNSRDC  ISSUES  TH3EE  TYPES  OF  REPORTS 

(1)  OTNSRDC  REPORTS.  A FORMAL  SERIES  PUBLISHING  INFORMATION  OF 
PERMANENT  TECHNICAL  VALUE,  DESIGNATED  BY  A SERIAL  REPORT  NUMBER. 

(21  DEPARTMENTAL  REPORTS.  A SEMIFORMAL  SERIES.  RECORDING  INFORMA- 
TION OF  A PRELIMINARY  OR  TEMPORARY  NATURE,  OR  OF  LIMITED  INTEREST  OR 
SIGNIFICANCE.  CARRYING  A DEPARTMENTAL  ALPHANUMERIC  1DENT*'  ICATIOM. 

(3)  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES,  USUALLY  '"TERN  _ 
WORKING  PAPERS  OR  DIRECT  REPORTS  TO  SPONSORS,  NUMBERED  AS  TM  SERIES 
REPORTS;  NOT  FOR  GENERAL  DISTRIBUTION. 


